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ABSTRACT 


The  theory  and  technique  of  Northrop  and  Thomsen  [1 980]  are 
applied  to  observations  of  energy  spectra  and  directional  aniso¬ 
tropies  of  0.6l  -  3.4l  MeV  protons  in  Saturn's  magnetosphere.  The 
observations  were  made  by  the  Goddard  Space  Flight  Center/University 
of  New  Hampshire  and  University  of  Iowa  instruments  aboard. 

Pioneer  11  during  the  Pioneer  encounter  with  Saturn  in  August/ 
September  1979-  Fourier  fits  to  15-minute  intervals  of  data  are 
combined  with  spectral  indices  to  yield  information  about  the  iT  x  f 
convection  velocity  and  temporal  changes  in  the  particle  population. 
There  is  a  fundamental  inability  to  distinguish  unambiguously 
between  the  two,  but  if  one  can  be  assumed,  the  other  then  follows 
from  these  calculations.  It  is  found  that  although  these  data  do 
not  by  themselves  allow  an  unambiguous  determination  of  the  extent 
of  corotation  in  Saturn's  outer  magnetosphere,  they  are  consistent 
with  exact  corotation  at  the  nominal  rotation  period  in  the 
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Introduction 

In  magnetospheric  physics,  the  extent  to  which  a  planet's 
ionosphere  is  electrically  coupled  to  its  magnetosphere  is  of  con¬ 
siderable  interest.  A  measure  of  this  coupling  is  the  degree  to 
which  the  charged  particles  in  the  magnetosphere  corotate  with  the 
planet.  One  thus  would  like  to  be  able  to  measure  the  2  X  5  con¬ 
vection  velocity  of  magnetospheric  particles  and  compare  it  with  the 
corotational  velocity. 

In  an  accompanying  paper,  Northrop  and  Thomsen  [1980]  outline 
the  theory  and  technique  by  which  spacecraft  observations  of  angular 
distributions  and  energy  spectra  of  energetic  particles  may  be  used 
to  deduce  2x3  flow  velocities  under  certain  circumstances.  In 
this  paper,  we  apply  the  results  of  Northrop  and  Thomsen  to  observa¬ 
tions  of  energetic  protons  in  Saturn's  magnetosphere  obtained  by  the 
University  of  Iowa  and  Goddard  Space  Flight  Center/University  of  New 
Hampshire  instruments  aboard  Pioneer  11. 

Figure  1,  from  Van  Allen  et  al.  [1980a]  shows  the  equatorial 
plane  projection  of  the  encounter  trajectory  of  Pioneer  11.  The 
plane  of  the  Pioneer  orbit  was  inclined  by  6.6°  relative  to  Saturn's 
equatorial  plane.  The  spacecraft  spin  axis  points  continuously  at 
the  Earth,  and  both  the  Iowa  and  GSFC/UNH  instruments  have  their 
view  direction  perpendicular  to  the  spacecraft  spin  axis  so  that  as 
the  spacecraft  spins,  the  detectors  scan  the  particle  distribution 
in  a  plane  perpendicular  to  the  spin  axis.  From  Figure  1  it  can  be 


seen  that  on  the  inbound  portion  of  the  encounter  trajectory,  the 
spacecraft  spin  axis  was  roughly  in  the  radial  direction  relative  to 
Saturn.  The  Iowa  and  GSFC/UNH  instruments  thus  scanned  in  a  plane 
nearly  perpendicular  to  the  radial  direction  and  were  therefore 
favorably  arranged  to  observe  anisotropies  associated  with  corota- 
tional  motion.  On  the  outbound  portion  of  the  trajectory,  the 
detectors  essentially  scanned  in  a  meridian  plane,  and  therefore 
azimuthal  flow  would  not  have  been  seen.  In  this  paper  we  report 
observations  made  on  the  inbound  portion  of  the  trajectory  between 
about  15  and  6  Rg. 
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Instrumentation 

The  University  of  Iowa  experiment  on  Pioneer  11  is  described 
elsewhere  in  this  issue  [Van  Allen  et  al.,  1980b].  In  this  paper  we 
are  concerned  with  observations  obtained  by  detector  G,  which  is  a 
thin  solid-state  detector  sensitive  to  protons  in  the  energy  range 
0.61  to  3-^1  MeV.  As  discussed  by  Van  Allen  et  al.  [1980b], 
detector  G  is  a  directional  detector  with  a  measured  angular 
response  function  having  FWHM  =  36°.  At  the  bit  rate  and  telemetry 
format  prevailing  during  the  Saturn  encounter,  detector  G  was 
sampled  twice  every  8.25  seconds,  each  sample  representing  an 
accumulation  period  of  0.75  seconds.  The  spacecraft  spin  period 
during  encounter  was  7-693  sec,  so  each  detector  readout  represented 
a  scan  through  35-1  degrees.  Each  sample  was  labelled  according  to 
the  spacecraft  roll  angle  of  the  detector  axis  at  the  center  of  the 
accumulation  period.  (The  roll  angle  is  the  angle  measured  in  the 
scan  plane  from  the  ascending  node  of  the  scan  plane  on  the 
ecliptic. ) 

Both  the  finite  angular  response  of  the  detector  and  the 
finite  sampling  time  have  the  effect  of  "smearing  out"  the  observed 
angular  distribution.  That  is,  if  the  true  angular  distribution  can 
oe  represented  by  a  Fourier  expansion  of  the  form 


R(cpl  =  M 


1  +  cos  i(cp  -  A^ ! 


i=l 


(1) 


6 


r 


then  a  Fourier  analysis  of  the  detector  counting  rate  would  yield 
smaller  values  for  the  Fourier  coefficients  than  the  true  FL  in  the 
actual  distribution.  This  reduction  factor  can  be  calculated  for 
each  Fourier  coefficient,  and  the  observed  values  can  thus  be  cor¬ 
rected  for  the  smear  effect.  The  correction  factors  for  the  Fourier 
coefficients  K^,  K^,  are  1.039,  1.160,  and  1.402,  Respectively. 

The  University  of  Iowa  detector  G  provides  good  angular 
resolution  but  no  spectral  information.  The  GSFC/UNH  instrument 
[Trainor  et  al.,  1974]  provides  detailed  energy  spectra  (0.2  <  < 

20  MeV)  spanning  the  energy  window  of  the  Iowa  detector.  The 
GSFC/UNH  experiment  also  measures  angular  distributions  which  agree 
closely  with  the  Iowa  observations.  Since  the  GSFC/UNH  experiment 
is  primarily  designed  for  cosmic  ray  observations,  its  telemetry 
assignment  during  encounter  was  minimal.  Therefore,  this  part  of 
the  GSFC/UNH  data  could  not  significantly  improve  the  overall 
statistical  accuracy  and  was  not  used.  The  two  experiments  thus 
provide  complementary  information,  and  we  combine  the  spectral 
information  available  from  the  Goddard  system  with  the  angular  dis¬ 
tributions  obtained  by  the  Iowa  detector  to  obtain  the  results  to  be 


presented  below. 
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Observations 

Counting  rates  from  the  University  of  Iowa  detector  G  have 
been  fit  in  15-minute  intervals  to  a  function  of  the  form  of  Equa¬ 
tion  (l),  with  i  =  1,2,3.  The  fitting  procedure  was  a  linear  least- 
squares  approach  rather  than  a  true  Fourier  analysis,  but  for  the 
nearly  uniform  roll  angle  coverage  represented  by  15  minutes  worth 
of  data,  the  two  procedures  are  equivalent.  Figure  2  shows  the 
results  of  the  fit  for  data  acquired  inbound  between  about  15  and 
5  Rg.  Note  that  there  is  a  significant  first  order  anisotropy  whose 
amplitude  (K^)  tends  generally  to  decrease  with  decreasing  distance 
from  Saturn  and  whose  phase  (A^  is  more  or  less  persistently 
parallel  to  the  equatorial  plane  (i.e.,  A^  ~  153°)  in  the  direction 
appropriate  for  corotation.  Such  a  first  order  anisotropy  is  quali¬ 
tatively  consistent  with  a  Compton- Get ting  effect  due  to  corotation. 
Note  also  the  strong  second,  order  anisotropy  whose  phase  (A2)  is 
also  approximately  parallel  to  the  equatorial  plane,  indicating  that 
these  are  "pancake"  distributions,  i.e.,  peaked  near  a  pitch  angle 
of  90°.  A  pancake  distribution  is  the  normal  signature  of  a  durably 
trapped  particle  population.  Figure  3  shows  sample  comparisons 
between  the  data  and  corresponding  fits  for  two  of  the  15-minute 
intervals . 

The  Fourier  coefficients  K^,  K^,  and  shown  in  Figure  2 
were  subsequently  corrected  for  the  smear  due  to  finite  detector 
geometry  and  sampling  time  as  discussed  above.  The  corrected  values 


are  related  to  the  Fourier  coefficients  (A^A^Ag, A^)  of  Northrop 
and  Thomsen  [1980]  by 

M  =  A 

o 

K  =  A. /A  i  =  1,2,3  .  (2) 

r  o 

The  phase  angles  Ai  (i  =  1,2,3)  shown  in  Figure  2  are  related  to  the 
phase  angles  x-  of  Northrop  and  Thomsen  [1980]  by 

h  =  Xi  +  (3) 

where  cpB  is  the  roll  angle  of  the  projection  of  the  magnetic  field 
vector  on  to  the  scan  plane  of  the  detector  (i.e.,  on  to  the  space¬ 
craft  equatorial  plane),  cp  is  thus  the  roll  angle  of  the  p  vector 

B 

in  the  notation  of  Northrop  and  Thomsen,  and  their  angle  \  is  just 
X  =  <P  -  <PB- 

With  the  notation  changes  discussed  above,  the  corrected 
Fourier  coefficients  and  phase  angles  can  be  directly  substituted 
into  Equations  (24)  and.  (25)  of  Northrop  and  Thomsen  [1980],  which 
relate  the  convection  velocity  Ug  to  the  spatial  gradient  of  the 
intensity  of  particles  having  an  energy  within  the  window  of  the 
detector.  For  ease  in  application,  Equations  (24)  and  (25)  can  be 
combined  to  show  the  dependence  on  each  Fourier  coefficient,  i.e., 
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2?  -  3 

7  -  1/2 


(y  +  l  -  k2) 


i 

l 


R7'1/2 


v  e  X  7M 

*  T  M—  -  KA  +  *2%  *  b% 
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where  y  is  the  differential  spectral  index  and 


s  .  ,-3L^g'(  i  -  »7;)x  s  *  .  s  )z-i 

'1  \  7  -  1  /  R7-3/2y[  1  sm  l  1  1 


1  -  ^  S,  x  ^  ^  ■' 


\  A  A 


2  -  cu  -1  -  9s  u>  1  &  +  \  PZ  -  |  cot  i  22/  •  ex 


T  .  *el  .  Sel  1  9el] 

X  [COt  1  9s  +  ei  X  9q  ■  sin  i  9p  J 

=»  (y  -  3/2V  1  -  R7”1  IT  *  X* 

^3  =  VT^i  Ax  .  r7-3/2/[“  n3  *  'qq  "  pp  +  cot  1  pzJ 


The  parameters  and  vectors  used  in  Equation  (4)  are  defined  in 
Northrop  and  Thomsen  [1980].  For  the  present  study,  the  unit  vector 
e^  in  the  direction  of  the  magnetic  field  (and  its  angle  of  inclina¬ 
tion  to  the  scan  plane,  i)  was  obtained  from  the  observed,  magnetic 
field  data  [courtesy  E.  J.  Smith  et  al.].  However,  the  components 
of  the  gradient  of  e^  (i.e.,  9e^/3s,  9e^/dp,  and  Se^/Bq)  were  cal¬ 
culated  assuming  a  dipole  field  model  [see  Equation  (23)  of  Northrop 
and  Thomsen].  The  dipole  approximation  is  applied  because  Smith 
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et  al.  [1980]  have  shown  that  during  the  Pioneer  11  encounter,  the 
dayside  magnetic  field  of  Saturn  was  indeed  quite  dipolar  in  shape, 
rather  than  being  flattened,  and  distended  by  a  large  current  sheet 
as  is  Jupiter's  magnetosphere.  For  the  University  of  Iowa  detec¬ 
tor  G,  the  ratio  R  of  lower  to  upper  threshold  energies  is  0.179? 
and  the  velocity  v^  of  protons  with  the  lower  threshold  energy  is 
1.08  X  10^  cm/sec.  The  local  proton  gyrofrequency  uu  was  calculated 
from  the  observed  magnetic  field  strength. 

The  final  observational  parameter  required  in  Equation  (4)  is 
the  differential  spectral  index  7.  Equation  (4)  is  based  on  the 
assumption  that  the  energy  spectrum  is  a  strict  power  law.  McDonald 
et  al.  [1980]  have  shown,  however,  that  beyond  about  11-5  Rg?  the 
proton  energy  spectrum  is  better  fit  by  a  sum  of  two  power  laws. 

The  left  panel  of  Figure  4  shows  an  example  of  such  a  spectrum.  We 
may  nevertheless  obtain  an  effective  v  from  such  a  spectrum  for  use 
in  Equation  (4)  because  the  fact  that  the  proton  intensity  decreases 
with  increasing  energy  means  that  an  integral  detector  is  princi¬ 
pally  responding  to  particles  with  energies  near  the  lower  detector 
threshold.  Thus  the  power  law  spectral  index  which  applies  in  the 
neighborhood  of  the  lower  threshold  energy  can  be  used  as  an  effec¬ 
tive  spectral  index  for  our  calculation.  Inside  of  11  Rc,  the 
spectrum  can  be  represented  by  a  single  power  law,  as  shown  in  the 
right  panel  of  Figure  4.  Figure  5  shows  the  spectral  index  deter¬ 
mined  by  the  GSFC/UNH  instrument  during  the  inbound  portion  of  the 
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encounter.  Each  point  is  located  at  the  center  of  the  30-minute 
interval  over  which  data  were  averaged  to  determine  the  spectrum. 
In  the  range  within  which  the  four-parameter  spectrum  provides  a 
better  fit  to  the  data,  the  effective  y  at  the  lower  threshold  of 
the  University  of  Iowa  detector  G  is  shown.  For  the  purpose  of 
evaluating  Equation  (4),  the  values  of  y  appropriate  to  the  time 
intervals  for  which  anisotropy  parameters  were  determined  were 
obtained  by  linearly  interpolating  between  the  values  plotted  in 
Figure  5- 
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Results 

In  Equation  (4)  it  can  be  seen  that  both  an  E?  X  $  (=  Ug) 
convection  velocity  and  a  spatial  gradient  of  the  density  (VM)  per¬ 
pendicular  to  the  magnetic  field  will  give  rise  to  a  first  order 
anisotropy.  For  a  spacecraft  near  the  magnetic  equatorial  plane 
such  as  Pioneer  11  at  Saturn,  a  radial  density  gradient  results  in 
an  ^  x  VM  contribution  to  the  anisotropy  which  is  in  the  azimuthal 
direction,  i.e.,  the  same  direction  as  the  anisotropy  arising  from  a 
corotational  Ug.  It  can  also  be  seen  from  Equation  (4)  that  a 
second  order  anisotropy,  in  the  presence  of  magnetic  curvature  or 
shear,  will  also  contribute  to  a  first  order  anisotropy.  There  will 
also  be  a  small  contribution  from  any  third  order  anisotropy,  but  as 
shown  in  Figure  2,  the  third  order  anisotropy  amplitude  is  generally 
quite  low. 

As  Northrop  and  Thomsen  [1980]  have  pointed  out,  the  separate 
contributions  of  spatial  gradients  and  E?  x  3  drifts  to  particle 
anisotropies  can  not  be  uniquely  identified  using  Equation  (4).  If, 
however,  one  knows  or  can  assume  something  about  the  spatial  gradi¬ 
ents,  Equation  (4)  allows  the  calculation  of  Ug.  Conversely,  if  Ug 
is  known,  the  gradient  can  be  calculated.  In  the  case  of  the 
Pioneer  11  encounter  with  Saturn,  there  are  two  different  situations 
we  may  explore:  One  is  a  steady-state  magnetosphere  and  the  other 
is  a  corotating  magnetosphere . 
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a.  Steady  State 

If  the  magnetospheric  configuration  and  particle  distribution 
were  unchanging  during  the  encounter,  then  the  changes  in  particle 
intensity  observed  by  Pioneer  11  as  it  traversed  the  magnetosphere 
(see  top  panel  of  Figure  2)  can  be  attributed  to  a  spatial  distribu¬ 
tion  alone.  Since  the  spacecraft  velocity  is  known,  the  gradient  of 
M  (at  least  the  component  of  the  gradient  parallel  to  the  spacecraft 
trajectory)  can  be  estimated.  Equation  (4)  then  yields  the  convec¬ 
tion  velocity  Ug. 

Figure  6  shows  the  azimuthal  component  of  Ug  which  results 
from  Equation  (4)  for  a  time- stationary  situation.  Since  the  space¬ 
craft  motion  was  predominantly  radial,  the  variations  in  M  seen  in 
Figure  2  have  been  attributed  to  a  purely  radial  gradient,  with  the 
nonradial  components  of  the  gradient  assumed  to  be  zero.  As  men¬ 
tioned  earlier,  a  radial  intensity  gradient  should  produce  an 
anisotropy  in  the  same  direction  as  a  corotational  u^.  The 
values  of  the  azimuthal  component  of  u^,  plotted  in  Figure  6 
have  been  smoothed  with  a  four-point  (i.e.,  one  hour)  running 
average  to  accentuate  lower-frequency  structure.  Also  shown  in 
Figure  6  as  a  solid  line  is  the  azimuthal  component  of  the  convec¬ 
tion  velocity  one  would  expect  for  particles  exactly  corotating 

ni 

with  the  nominal  planetary  rotation  period  at  10  4o  [Kaiser 
et  al. ,  I960]. 


A  weighted  least-squares  fit  of  the  azimuthal  velocity  com¬ 


ponent  to  a  function  of  the  form  u^,  =  rH  gives 


Q  =  [1.729  (±  0.229)]  X  10~k  s'1 


which  corresponds  to  a  rotational  period  of 


P 
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This  rotation  period  gives  the  dashed  line  in  Figure  6.  Thus, 
although  Figure  6  shows  apparently  significant  deviations  from  exact 
corotation  when  a  steady  state  is  assumed,  on  the  average  the 
results  are  consistent  with  corotation  at  the  nominal  rate. 


b.  Exact  Corotation 

If  the  magnetosphere  is  corotating  exactly,  then  u^,  is 
known  and  Equation  (4)  yields  the  gradient  of  the  particle 
intensity.  This  gradient,  dotted  with  the  spacecraft  velocity, 
tells  us  the  convective  part  of  the  time  variations  in  the 
intensity  observed  by  the  spacecraft.  Subtracting  this  from 
the  total  time  derivative  actually  observed  by  the  spacecraft 
leaves  an  estimate  of  the  purely  temporal  changes  in  the 
intensity,  i.e.,  d(log  M)/dt. 
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Figure  7  shows  the  values  of  ^(log  M)/dt  which  result  from 
assuming  exact  corotation  at  the  nominal  period  of  10h  4om  (again 
smoothed  with  a  four-point  running  average).  There  are  intervals  of 
time  during  which  the  partial  time  derivative  is  found  to  be  sig¬ 
nificantly  nonzero,  but  a  weighted  average  over  all  the  data  shown 
in  Figure  7  gives  an  average  time  change  over  the  entire  interval  of 

{d'Uft' =  [_  4.9  (±  i4.5)]  x  io'U  min-1 

The  average  temporal  change  when  exact  corotation  is  assumed  is  thus 
zero  within  the  errors. 

To  test  the  sensitivity  of  this  result,  we  have  also  computed 
^(log  M)/dt  for  the  case  Ug  =  0.  These  results  are  shown  in 
Figure  8.  As  can  be  seen  from  this  figure,  if  zero  corotation  is 
assumed,  the  deduced  partial  time  derivative  is  negative  at  almost 
all  times.  The  average  over  the  entire  interval  is 

{^ft  =  t'  1,07  0-21)]  X  ID'2  min-1  , 

which  would  represent  a  net  decrease  in  the  energetic  particle 
intensity  by  roughly  a  factor  of  650  during  this  time  period. 

Such  a  large  decrease  in  the  energetic  particle  content  of 
the  outer  magnetosphere  seems  too  large  to  be  credible. 
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Discussion 

As  discussed  in  Northrop  and  Thomsen  [1980],  Equation  (U) 
is  derived  on  the  basis  of  a  model  gyrotropic  distribution  function 
(f  in  their  notation)  which  is  pancake  in  shape,  i.e. ,  peaked 
near  a  pitch  angle  of  90°.  This  model  is  appropriate  for  the 
present  study  because,  as  shown  in  Figure  2,  the  energetic  proton 
distribution  was  persistently  pancake  in  nature  (a2  ~  153°) 
throughout  Saturn's  dayside  magnetosphere.  The  model  distribution 
function  f  also  did  not  include  a  provision  for  parallel  flow, 
but  as  can  be  seen  in  Figure  2,  the  exclusion  of  parallel  flow 
from  the  treatment  of  anisotropies  at  Saturn  appears  to  be  valid 
since  the  first  order  anisotropy  direction  A-^  is  always  observed 
to  be  quite  close  to  a  direction  perpendicular  to  the  magnetic 
field.  (A  field-aligned  flow  would  have  A-^  ~  63°  or  2U5°.) 

This  is  in  marked  contrast  to  the  situation  during  the  Pioneer  10 
encounter  with  Jupiter,  where  Northrop  et  al.  [1979]  found  that 
the  effects  of  parallel  flow  were  very  important  and  could  not 
be  excluded. 

Figure  6  shows  that  on  the  average  the  convection  velocity 
derived  from  Equation  (U)  when  temporal  changes  are  ignored  is 
consistent  with  corotation  at  the  nominal  period.  However,  there 
are  also  times  during  which  the  apparent  velocity  would  be  sig¬ 
nificantly  different  from  the  corotation  value.  Of  the  38  points 
shown  in  Figure  6,  22  (or  58 %)  of  them  deviate  from  the  value 


expected  for  exact  corotation  by  more  than  la,  3  (or  8%)  lie  at 
least  2 a  away  from  exact  corotation,  and  1  (or  3$)  lies  at  least 
3a  away.  For  a  normal  distribution  of  errors,  only  32 5%,  and 
<  1^  should  lie  more  than  la,  2 a,  and  3a,  respectively,  away  from 
the  mean.  A  situation  such  as  that  near  11  Rg,  where  the  deduced 
velocity  is  more  than  twice  the  value  expected  for  exact  corotation, 
would  appear  to  be  physically  unrealistic,  so  the  assumption  of 
steady  state  is  probably  inappropriate  during  at  least  part  of 
this  period. 

Large-scale  temporal  changes  are  not  unreasonable  in  light  of 
the  highly  disturbed  solar  wind  conditions  that  were  present  during 
the  Pioneer  11  Saturn  encounter.  Wolfe  et  al.  [1980]  have  reported 
that  the  upstream  solar  wind,  dynamic  pressure  was  as  much  as  7  to  8 
times  as  high  as  its  normal  value  during  the  period  just  before  the 
inbound  encounter  with  the  magnetosphere.  They  attribute  the  second 
inbound  bow  shock  crossing  to  a  temporary  increase  in  solar  wind 
dynamic  pressure,  which  apparently  pushed  the  bow  shock  back  past 
Pioneer  11,  leaving  it  again  in  the  solar  wind  until  it  re¬ 
encountered  the  bow  shock  some  5  hours  later,  more  than  1%  closer 
to  the  planet.  By  the  time  Pioneer  11  exited  from  the  magnetosphere 
on  the  dawnside  3  l/2  days  later,  the  solar  wind  dynamic  pressure 
had  diminished  to  near  or  slightly  below  normal.  Thus,  there  is 
good  evidence  that  the  size  of  Saturn's  magnetosphere  was  probably 
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varying  considerably  during  the  Pioneer  11  encounter.  Such  varia¬ 
tions  in  the  magnetospheric  field  are  known  to  result  in  redistribu¬ 
tion  and  energization  of  trapped  particles  [see,  e.g.,  Brewer  et 
al.,  1969]-  Thus,  we  should  not  be  surprised  that  the  steady-state 
assumption  results  in  improbable  convection  velocities. 

Granting  that  temporal  changes  were  occurring  during  the 
period  under  consideration,  we  no  longer  can  estimate  point-by-point 
the  I  x  I  velocity  on  the  basis  of  Pioneer  11  observations.  We  can, 
however,  ask  whether  or  not  the  assumption  of  corotation  leads  to 
unreasonable  values  for  the  temporal  changes  that  are  probably 
occurring.  Referring  to  Figure  7>  we  see  that  there  are  several 
intervals  of  time  during  which  the  derived  partial  time  derivative 
of  the  intensity  is  persistently  either  positive  or  negative. 
Integrating  d(log  M)/dt  over  such  an  interval  tells  us  the  total 
temporal  change  in  M  during  that  period.  For  example,  the  interval 
of  persistently  negative  partial  time  derivative  from  05*+  5  to  0730 
represents  a  net  decrease  in  the  intensity  by  a  factor  of  only  about 
0.46,  and  the  interval  of  persistently  positive  partial  time  deriva¬ 
tive  from  1000  to  1115  represents  a  net  increase  in  M  by  a  factor  of 
1.7*  Thus,  the  temporal  changes  deduced  from  the  assumption  of 
exact  corotation  are  always  less  than  a  factor  of  2  one  way  or  the 
other,  and  as  we  noted  earlier,  the  net  change  over  the  entire 
period  we  have  studied  here  is  zero.  On  the  other  hand,  the  assump¬ 
tion  of  zero  corotation  leads  to  persistently  negative  time 
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derivatives  throughout  the  entire  interval  and  a  net  decrease  in  the 
overall  intensity  by  almost  three  orders  of  magnitude,  Such  a  large 
change  seems  highly  improbable  even  in  light  of  the  very  disturbed 
interplanetary  conditions. 

The  error  bars  shown  on  the  data  points  in  Figures  6-8  have 
been  computed  by  propagating  the  errors  in  the  and  A^  through 
Equation  (4),  assuming  the  errors  in  the  and  A^  are  independent. 
We  have  also  tested  the  sensitivity  of  our  results  to  the  adopted 
value  of  the  spectral  index  by  performing  the  calculations  for  a 
spectral  index  equal  to  the  ..ominal  value  ±  0.5.  Such  a  variation 
in  7  results  in  derived  azimuthal  velocities  which  vary  by  less 
than  about  15$  from  those  shown  in  Figure  6  and  which  lie  within 
the  error  bars  shown  in  that  figure.  Thus,  although  uncertainties 
if  that  magnitude  in  the  spectral  index  might  slightly  affect  our 
quantitative  results,  our  basic  conclusions  would  not  be  signifi¬ 
cantly  affected. 

Iri  summary,  although  these  data  by  themselves  do  not  allow 
an  unambiguous  determination  of  the  extent  of  corotation  in  Saturn  s 
outer  magnetosphere,  they  are  consistent  with  exact  corotation  at 
the  nominal  rotation  period  in  the  presence  of  significant  but  not 
unreasonable  temporal  variations  in  the  proton  population.  Such  a 
picture  is  supported  by  data  from  the  plasma  analyzer  on  Fioneer  11 
which,  as  reported  by  Wolfe  et  al .  [1980]  and  Frank  et  al.  [this 
issue!,  found  plasma  flow  in  Saturn's  magnetosphere  which  appears  to 
>,0  -nuni stent  with  corotation. 
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FIGURE  CAPTIONS 


1.  Equatorial  plane  projection  of  the  encounter  trajectory 
of  Pioneer  11  at  Saturn.  The  spacecraft  spin  axis  points 
continuously  at  Earth,  and  the  detectors'  look  direction  is 
perpendicular  to  the  spin  axis. 

2.  Fourier  coefficients  and  phase  angles  from  fits  to 
15-minute  intervals  of  data  from  the  University  of  Iowa 
detector  G  during  the  inbound  portion  of  the  Pioneer 
encounter  with  Saturn.  Error  bars  axe  calculated  from  the 
RMS  deviation  between  the  fit  and  the  data. 

5.  Comparisons  between  15-minute  intervals  of  data  and  the 
corresponding  Fourier  fits  as  described  in  the  text.  Repre¬ 
sentative  error  bars  show  the  statistical  uncertainty  in  the 
individual  points. 

h ■  Sample  fits  to  the  proton  energy  spectrum.  (Left) 


Adopted  spectral  form  is  the  sum  of  two  power  laws.  (Right) 
Adopted  spectral  form  is  a  single  power  law. 
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Figure  5-  Effective  power  law  differential  spectral  index  at 

Ep  =  0.6l  MeV.  Values  at  times  between  plotted  points  are 
obtained  by  linear  interpolation. 

Figure  6.  Azimuthal  component  of  the  E?  x  ^  convection  velocity 
deduced  from  Equation  (4)  for  a  time-stationary  magneto- 
spheric  configuration.  The  points  have  been  smoothed  by  a 
running  four-point  average.  The  error  bars  are  thus  approxi¬ 
mately  half  the  size  of  the  estimated  errors  in  the  original 
unsmoothed  15-minute  results.  Errors  have  been  estimated  by 
propagating  the  errors  in  the  and  through  Equation  (4), 
assuming  the  errors  in  the  and  are  independent.  Also 

shown  is  the  azimuthal  component  of  the  convection  velocity 

h  m 

one  would  expect  for  corotation  with  a  period  of  10  40 

(solid  line)  and  the  least-squares  fit  to  the  data  points 
(dashed  line),  which  corresponds  to  a  period  of  10*1  6m. 

Figure  J.  Partial  time  derivative  of  the  logarithm  of  the  proton 

intensity  deduced  from  Equation  (4)  for  a  magnetosphere  which 
corotates  with  a  period  of  10^  4om.  The  points  have  been 
smoothed  by  a  running  four-point  average,  and.  the  comments  in 
Figure  6  regarding  the  error  estimates  are  applicable  here  as 
well . 

Figure  tf.  fame  as  Figure  7  except  the  magnetosphere  is  assumed  not 
* o  be  corotating  at  all,  i.e.,  u  =  0. 
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